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Abstract. Nonclassical states of light play a central role in many quan¬ 
tum information protocols. Very recently, their quantum features have 
been exploited to improve the readout of information from digital memo¬ 
ries, modeled as arrays of microscopic beam splitters [S. Pirandola, Phys. 
Rev. Lett. 106, 090504 (2011)]. In this model of “quantum reading”, a 
nonclassical source of light with Einstein-Podolski-Rosen correlations has 
been proven to retrieve more information than any classical source. In 
particular, the quantum-classical comparison has been performed un¬ 
der a global energy constraint, i.e., by fixing the mean total number of 
photons irradiated over each memory cell. In this paper we provide an 
alternative analysis which is based on a local energy constraint, mean¬ 
ing that we fix the mean number of photons per signal mode irradiated 
over the memory cell. Under this assumption, we investigate the critical 
number of signal modes after which a nonclassical source of light is able 
to beat any classical source irradiating the same number of signals. 


1 Introduction 

Quantum information has disclosed a modern approach to both quantum me¬ 
chanics and information theory [1]. Very recently, this field has been developed 
into the so-called “continuous variable” setting, where information is encoded 
and processed using quantum systems with infinite dimensional Hilbert spaces 
[2-6]. Bosonic systems, such as the radiation modes of the electromagnetic field, 
are today the most studied continuous variable systems, thanks to their strict 
connection with quantum optics. In the continuous variable framework, a wide 
range of results have been successfully achieved, including quantum telepor¬ 
tation [7-12], teleportation networks [13-15] and games [16,17], entanglement 
swapping protocols [18-20], quantum key distribution [21-25], two-way quantum 
cryptography [26,27], quantum computation [28-34] and cluster quantum com¬ 
putation [35-38]. Other studies have lead to the full classification of Gaussian 
channels and collective Gaussian attacks [39-41], the computation of secret-key 
capacities and their reverse counterpart [42-46], and possible schemes for quan¬ 
tum direct communication [47,48]. 

One of the key resources in many protocols of quantum information is quan¬ 
tum entanglement. In the bosonic setting, quantum entanglement is usually 
present under the form of Einstein-Podolski-Rosen (EPR) correlations [49], where 


the quadrature operators of two separate bosonic modes are so correlated to beat 
the standard quantum limit [50]. The simplest source of EPR correlations is the 
two-mode squeezed vacuum (TMSV) state. In the number-ket representation 
this state is defined by 

OO 

10 = (cosh^-^ ^(tanh^” |n)Jn)^ , 

n—0 

where ^ is the squeezing parameter and {s,i} is an arbitrary pair of bosonic 
modes, that we may call “signal” and “idler”. In particular, ^ quantifies the 
signal-idler entanglement and gives the mean number of photons sinh^^ in each 
mode. Since it is entangled, the TMSV state cannot be prepared by applying 
local operations and classical communications (LOCCs) to a couple of vacua 
|0)s ® O’" other kind of tensor product state. For this reason, the TMSV 

state cannot be expressed as a classical mixture of coherent states |/3)j with 
a and /3 arbitrary complex amplitudes. In other words, its P-representation [58, 
59] 

\0{^\=J j (facfp'P{a,l3) |a)s (a| 0 |/3)i (/3| , 

involves a function V which is non-positive and, therefore, cannot be consid¬ 
ered as a genuine probability distribution. For this reason, the TMSV state is 
a particular kind of “nonclassical” state. Other kinds are single-mode squeezed 
states and Fock states. By contrast a bosonic state is called “classical” when its 
P-representation is positive, meaning that the state can be written as a classical 
mixture of coherent states. Thus a classical source of light is composed by a set 
of m bosonic modes in a state 

P = J (famViai,--- ,am) 0™^! ]«*) («*] , (I) 

where V is positive and normalized to I. Typically, classical sources are just 
made by collection of coherent states with amplitudes {di,--- ,dm}, be., p = 
®'k=i \^k) (dfej which corresponds to have 

m 

^ S‘^{ak - ak) . 
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In other situations, where the sources are particularly chaotic, they are better de¬ 
scribed by a collection of thermal states with mean photon numbers {ri-i, • • • , rim}, 
so that 

exp(-lafcl rife) 

TThfc 

More generally, we can have classical states which are not just tensor products 
but they have (classical) correlations among different bosonic modes. 

The comparison between classical and nonclassical states has clearly triggered 
a lot of interest. The main idea is to compare the use of a candidate nonclassical 


m 

p=n 



state, like the EPR state, with all the classical states for specific information 
tasks. One of these tasks has been the detection of low-reflectivity objects in far 
target regions under the condition of extremely low signal-to-noise ratios. This 
scenario has been called “quantum illumination” and has been investigated in a 
series of papers [51-56]. 

Most recently, the EPR correlations have been exploited for a completely 
different task in a completely different regime of parameters. In the model of 
“quantum reading” [57], the EPR correlations have been used to retrieve in¬ 
formation from digital memories which are reminiscent of today’s optical disks, 
such as CDs and DVDs. A digital memory can in fact be modelled as a sequence 
of cells corresponding to beam splitters with two possible reflectivities ro and 
ri (used to encode a bit of information). By fixing the mean total number of 
photons N irradiated over each memory cell, it is possible to show that a non- 
classical source of light with EPR correlations retrieves more information than 
any classical source [57]. In general, the improvement is found in the regime of 
few photons {N = I -y 100) and for memories with high reflectivities, as typi¬ 
cal for optical memories. In this regime, the gain of information given by the 
quantum reading can be dramatic, i.e., close to I bit for each bit of the memory. 

An important point in the study of Ref. [57] is that the quantum-classical 
comparison is performed under a global energy constraint, i.e., by fixing the total 
number of photons N which are irradiated over each memory cell (on average). 
Under this assumption, it is possible to construct an EPR transmitter, made by a 
suitable number of TMSV states, which is able to outperform any classical source 
composed by any number of modes. In the following we consider a different and 
easier comparison: we fix the number of signal modes irradiated over the target 
cell (M) and the mean number of photons per signal mode (Ns). Under these 
assumptions, we compare an EPR transmitter with a classical source. Then, for 
fixed Ns, we determine the critical number of signal modes after which 

an EPR transmitter with M > is able to beat any classical source (with 

the same number of signals M). 

2 Readout mechanism 

Here we briefly review the basic readout mechanism of Ref. [57], specifying the 
study to the case of a local energy constraint. Let us consider a model of digital 
optical memory (or disk) where the memory cells are beam splitter mirrors with 
different reflectivities r = ro, ri (with ri > rp). In particular, the bit-value u = 0 
is encoded in a lower-reflectivity mirror (r = ro), that we may call a pit, while 
the bit-value u = I is encoded in a higher-reflectivity mirror (r = ri), that 
we may call a land (see Fig. I). Close to the disk, a reader aims to retrieve 
the value of the bit u which is stored in each memory cell. For this sake, the 
reader exploits a transmitter (to probe a target cell) and a receiver (to measure 
the corresponding output). In general, the transmitter consists of two quantum 
systems, called signal S and idler I, respectively. The signal system A is a set of 
M bosonic modes which are directly shined on the target cell. The mean total 



number of photons of this system is simply given by fV = MNs, where Ns is the 
mean number of photons per signal mode (simply called “energy”, hereinbelow). 
At the output of the cell, the reflected system R is combined with the idler 
system /, which is a supplementary set of bosonic modes whose number L can 
be completely arbitrary. Both the systems R and I are finally measured by the 
receiver (see Fig. 1). 


Thermal 



Fig. 1. Model of memory. Digital information is stored in a disk whose memory 
cells are beam splitter mirrors with different reflectivities: r — ro encoding bit-value 
M = 0 and r = n encoding bit-value u — 1. Readout. A reader is generally composed 
by a transmitter and a receiver. It retrieves a stored bit by probing a memory cell with 
a signal system S (M bosonic modes) and detecting the reflected system R together 
with an idler system I {L bosonic modes). In general, the output system R combines 
the signal system S with a bath system B (M bosonic modes in thermal states). The 
transmitter is in a state p which can be classical (classical transmitter) or non-classical 
(quantum transmitter). In our work, we consider a quantum transmitter with EPR 
correlations between signal and idler systems. 


We assume that Alice’s apparatus is very close to the disk, so that no signifi¬ 
cant source of noise is present in the gap between the disk and the decoder. How¬ 
ever, we assume that non-negligible noise comes from the thermal bath present 
at the other side of the disk. This bath generally describes stray photons, trans¬ 
mitted by previous cells and bouncing back to hit the next ones. For this reason, 
the reflected system R combines the signal system S with a bath system B oi M 
modes. These environmental modes are assumed in a tensor product of thermal 
states, each one with Nb mean photons (white thermal noise). In this model 
we identify five basic parameters: the reflectivities of the memory {ro,ri}, the 
temperature of the bath Nb, and the profile of the signal {M,Ns}, which is 
given by the number of signals M and the energy Ns- 

In general, for a fixed input state p at the transmitter (systems S,I), Alice 
will get two possible output states ctq and <Ji at the receiver (systems R,I). 
These output states are the effect of two different quantum channels, fo and £i, 
which depend on the bit u = 0,1 stored in the target cell. In particular, we have 
cr„ = {£u ®X){p), where the conditional channel £u acts on the signal system, 
while the identity channel I acts on the idler system. More exactly, we have 
where is a one-mode lossy channel with conditional loss and 
fixed thermal noise Nb- Now, the minimum error probability Perr affecting the 



decoding of u is just the error probability affecting the statistical discrimination 
of the two output states, cto and cti, via an optimal receiver. This quantity is 
equal to Perr = [1 ~ CTi)]/2, where D(ao, cri) is the trace distance between 

(To and CTi [60-62]. Clearly, the value of Perr determines the average amount of 
information which is decoded for each bit stored in the memory. This quantity 
is equal to J = 1 — H{Perr), where H{x) := —a;log 2 x — {1 — x) log 2 (l — x) is 
the usual formula for the binary Shannon entropy. In the following, we com¬ 
pare the performance of decoding in two paradigmatic situations, one where the 
transmitter is described by a non-classical state (quantum transmitter) and one 
where the transmitter is in a classical state (classical transmitter). In particular, 
we show how a quantum transmitter with EPR correlations (EPR transmitter) 
is able to outperform classical transmitters. The quantum-classical comparison 
is performed for fixed signal profile {M,Ns}. Then, for various fixed values of 
the energy Ns (local energy constraint), we study the critical number of signal 
modes after which an EPR transmitter (with M > signals) is able 

to beat any classical transmitter (with the same number of signals M). 


3 Quantum-classical comparison 


First let us consider a classical transmitter. A classical transmitter with M 
signals and L idlers is described by a classical state p as specified by Eq. (I) with 
m = M L.lxi other words it is a probabilistic mixture of multi-mode coherent 
states \cik) {oik\- Given this transmitter, we consider the corresponding 

error probability Pelr^^ which affects the readout of the memory. Remarkably, 
this error probability is lower-bounded by a quantity which depends on the 
signal profile {M, Ns}, but not from the number L of the idlers and the explicit 
expression of the P-function. In fact, we have [57] 


P:^r > C{M, Ns) := 1 ^ ( 2 ) 

where F{Ns) is the fidelity between 77.o(|A^y^)(iVy^|) and TZi{\Ny^){N^^ \), 
the two possible outputs of the single-mode coherent state \Ny^){Ny^\. As a 
consequence, all the classical transmitters with signal profile {M,Ns} retrieve 
an information which is upper-bounded by Jdass := I — H[C{M, Ns)]. 

Now, let us construct a transmitter having the same signal profile {M,Ns}, 
but possessing EPR correlations between signals and idlers. This is realized by 
taking M identical copies of a TMSV state, i.e., p = |^) where Ns = 

sinh^^. Given this transmitter, we consider the corresponding error probability 
affecting the readout of the memory. This quantity is upper-bounded by 
the quantum Chernoff bound [63-67] 


< Q{M,Ns) 
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where := (77.„ (8)I)(|^)(^|). Since 9q and 9i are Gaussian states, we can 
write their symplectic decompositions [68] and compute the quantum Chernoff 



bound using the general formula for multimode Gaussian states given in Ref. [67]. 
Then, we can easily compute a lower bound Jquant := 1 — H[Q{M,Ns)] for the 
information which is decoded via this quantum transmitter. 

In order to show an improvement with respect to the classical case, it is suffi¬ 
cient to prove the positivity of the “information gain” G := Jquant — Jclass- This 
quantity is in fact a lower bound for the average information which is gained 
by using the EPR quantum transmitter instead of every classical transmitter. 
Roughly speaking, the value of G estimates the minimum information which is 
gained by the quantum readout for each bit of the memory. In general, G is a 
function of all the basic parameters of the model, i.e., G = G(M, iVg, ro, ri, TVs). 
Numerically, we can easily find signal profiles {M, Ns}, classical memories {rp, ri}, 
and thermal baths Nb, for which we have the quantum effect G > 0. Some of 
these values are reported in the following table. 


M 

Ns 

ro 

ri 

Nb 

G (bits) 

1 

3.5 

0.5 

0.95 

0.01 

6.2 X 10"^ 

10 

1 

0.2 

0.8 

0.01 

3.4 X 10"^ 

30 

1 

0.38 

0.85 

1 

1.2 X 10"^ 

100 

0.1 

0.25 

0.85 

0.01 

5.9 X 10"^ 

200 

0.1 

0.6 

0.95 

0.01 

0.22 


0.01 

0.995 

1 

0 



Note that we can find choices of parameters where G ~ 1, i.e., the classical 
readout of the memory does not decode any information whereas the quantum 
readout is able to retrieve all of it. As shown in the last row of the table, this 
situation can occur when the both the reflectivities of the memory are very close 
to I. From the first row of the table, we can acknowledge another remarkable 
fact: for a land-reflectivity ri sufficiently close to I, one signal with few photons 
can give a positive gain. In other words, the use of a single, but sufficiently 
entangled, TMSV state |^) (^| can outperform every classical transmitter, which 
uses one signal mode with the same energy (and potentially infinite idler modes). 

According to our numerical investigation, the quantum readout is generally 
more powerful when the land-reflectivity is sufficiently high (i.e., ri > 0.8). For 
this reason, it is very important to analyze the scenario in the limit of ideal 
land-reflectivity (ri = I). Let us call “ideal memory” a classical memory with 
ri = I. Clearly, this memory is completely characterized by the value of its pit- 
reflectivity ro- For ideal memories, the quantum Chernoff bound of Fq. (3) takes 
the analytical form 

Q = ^{[1 + (1 - V^)Nsf + NBi2Ns + 1)(1 - ro)}-^, 
and the classical bound of Eq. (2) can be computed using 
F{Ns) = 7 "^exp[- 7 "^(I - y^fNs] , 

where 7 := 1 -I- (I — ro)NB [57]. Using these formulas, we can study the behavior 
of the gain G in terms of the remaining parameters {M, Ns,ro, Nb}- Let us 











































consider an ideal memory with generic ro G [0,1) in a generic thermal bath Nb > 
0. For a fixed energy Ns, we consider the minimum number of signals 
above which G > 0 [69]. This critical number can be defined independently from 
the thermal noise Nb (via an implicit maximization over Nb)- Then, for a given 
value of the energy Ns, the critical number is a function of tq alone, i.e., 

][f(Ns) — Its behavior is shown in Fig. 2 for different values of the 

energy. 



Fig. 2. Number of signals M (logarithmic scale) versus pit-reflectivity ro. The curves 
refer to Ns = 0.01, 0.1 and 0.5 photons. For each value of the energy Ns, we plot 
the critical number as function of ro. All the curves have an asymptote 

at ro = 1. For Ns > 2.5 photons (curves not shown), we have another asymptote at 
ro = 0. 


It is remarkable that, for low-energy signals {Ns = O.OI G- I photons), the 
critical number M^^^'^{ro) is finite for every rp G [0,1). This means that, for 
ideal memories and low-energy signals, there always exists a finite number of 
signals above which the quantum readout of the memory is more efficient 

than its classical readout. In other words, there is an EPR transmitter with 
M > able to beat any classical transmitter with the same number of 

signals M. In the considered low-energy regime, M(^®)(ro) is relatively small 
for almost all the values of tq, except for tq —f 1 where M(^®)(ro) — f oo. In fact, 
for To ~ 1, we derive M(^®)(ro) — [ 4 A^ 5 ( 2 A ^5 -h I)(I — ro)]“^, which diverges at 
ro = 1. Such a divergence is expected, since we must have = 1/2 

for ro = ri (see Appendix A for details). Apart from the divergence at ro = 1, 
in all the other points vq G [0,1), the critical number M(^®^(ro) decreases for 
increasing energy Ns (see Fig. 2). In particular, for Ns = 1 photon, we have 
M('^®)(ro) ~ 1 for most of the reflectivities tq- In other words, for energies 
around one photon, a single TMSV state is sufficient to provide a positive gain 
for most of the ideal memories. However, the decreasing trend of M(^®^(ro) 
does not continue for higher energies {Ns > 1). In fact, just after Ns = 1, 
M(^®)(ro) starts to increase around tq = 0. In particular, for Ns > 1, we can 
derive M(^®^( 0 ) ~ (ln 2 )[ 21 n(l -|- Ns) — which is increasing in Ns, and 

becomes infinite at Ns — 2.5. As a consequence, for Ns > 2.5 photons, we have a 
second asymptote appearing at ro = 0 (see Appendix B for details). This means 




that the use of high-energy signals {Ns > 2.5) does not assure positive gains for 
memories with extremal reflectivities ro = 0 and ri = 1. 


4 Conclusion 

In conclusion, we have considered the basic model of digital memory studied in 
Ref. [57], which is composed of beam splitter mirrors with different reflectivi¬ 
ties. Adopting this model, we have compared an EPR transmitter with classical 
sources for fixed signal profiles, finding positive information gains for memo¬ 
ries with high land-reflectivities (ri > 0.8). Analytical results can be derived 
in the limit of ideal land-reflectivity (ri = 1) which defines the regime of ideal 
memories. In this case, by fixing the mean number of photons per signal mode 
(local energy constraint), we have computed the critical number of signals after 
which an EPR transmitter gives positive information gains. For low-energy sig¬ 
nals (0.01-^ I photons) this critical number is finite and relatively small for every 
ideal memory. In particular, an EPR transmitter with one TMSV state can be 
sufficient to achieve positive information gains for almost all the ideal memories. 
Finally, our results corroborate the outcomes of Ref. [57] providing an alternative 
study which considers a local energy constraint instead of a global one. As dis¬ 
cussed in Ref. [57], potential applications are in the technology of optical digital 
memories and involve increasing their data-transfer rates and storage capacities. 


A General asymptote at ro = 1 


According to Fig. 2, the critical number M^^^'>{ro) diverges for tq —>■ 1. Let us 
analyze the behavior of G around the singular point tq = 1, by setting ro = 1 — £ 
and expanding G for £ —>• 0+. It is easy to check that, for every we have 
G > 0 if and only if M > [4A^5(2A^5 -I- 1)£]“^. In particular, in the absence of 
thermal noise {Nb =0), we have 


MNs{AMNs - 1)£2 , 

o 1 ^ o 


(4) 


which is positive if and only if M > {4:Ns)~^. 

Proof. Note that G > 0 if and only if A := Q{M, Ns) —C{M, Ns) < 0. Let 
us expand A = A{M, Ns, Nb, 1 — £) at the first order in £. For a given Nb > 0, 


(MNbs)^^^ 
negative if and only if 


we have A = ^ 


M{Nb +Ns + 2NBNs)e 


+ 0(£^/^), which is 


M > 


Nb 

{Nb + Ns + 2 NbNs)'^s 


k{Nb) ■ 


Notice that k{Nb) is maximum for Ng = Ns{l + 2Ns) Then, for every 
Nb > 0, we have Z\ < 0 if and only if 

I 


M > k{N^) 


ms{2Ns + 1 )£ ■ 


(5) 



Now, let us consider the particular case of Nb = 0. In this case, we have the 
hrst-order expansion A = [1~2 (MA^s)^^^]e/4+0(£^), or equivalently 

the second-order expansion of G given in Eq. (4). It is clear that Z\ < 0, i.e., 
G > 0, when M > I/4iV5. However, this condition is less restrictive than the 
one in Eq. (5) which, therefore, can be extended to every Nb > 0. ■ 


B High-energy asymptote at rq = 0 


Let us analyze the behavior of for Ns > I and vq = 0. One can check 

that, for Ns > 1, the greatest value of occurs when Nb = 0. In this 

case, i.e., for tq = Nb = 0 and ri = 1, we have Q{M, Ns) = [(1 -|- Ns)~‘^^]/2, 
and 


C{M,Ns) 


1 _ VI - 


:=C° 


Let us consider the critical value of M such that G{M, Ns) = 0, which 

is equivalent to Q = C. We also consider the value M such that Q = C°°. Since 
C°° < C, we have that > M. Actually, we find that ~ M 

with very good approximation when Ns > 1 (see Eig. 3). Then, for every Ns > 1, 
we can set (0) ~ (In 2) [2 ln(I -|- Ns) — The latter quantity becomes 

infinite for 21n(I -|- Ns) = Ns, i.e., for Ns > 2.51 photons. 



Ns 

Fig. 3. Minimum number of signals M versus energy Ns- The solid curve represents 
M{^®)(0) while the dashed curve represents M. Note that the minimum number of 
signals is actually given by [M] where [• • • ] is the ceiling function. 
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